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ABSTRACT 


The line positions, line intensities, and collisional line half-widths have 
been calculated for the fundamental, first-overtone, and second- overtone 
vibration-rotation bands of C 12 0 16 and C 13 0 16 . The C 12 0 16 and C 13 O ie line 
intensities were computed for temperatures ranging from 175 to 3500°K and 
were weighted in accordance with a terrestrial C 12 /C 13 abundance ratio. All 
lines within 5 to 6 orders of magnitude of the maximum line intensity of a 
band have been considered (~4500 lines at 3500°K). Tables of the above mo- 
lecular absorption line parameters are being published (Kunde, 1967). 

The CO mass absorption coefficient has been computed for several late 
spectral type stellar -atmosphere conditions. The straight mean absorption 
coefficient has been evaluated for these conditions in several ways: (1) ex- 
actly, by using line-by-line integration techniques; (2) approximately, by using 
the quasi-random molecular-absorption band model; and (3) approximately, 
by computing the expression S T m /Av, where S T m is the total line intensity in 
spectral interval Av. Results of this investigation indicate that the simpler 
S T m /Ai> approximation represents the CO mean absorption coefficient in the 
late type atmospheres generally to an accuracy better than 1 percent, if the 
spectral intervals are greater than 5 cm' 1 . The quoted accuracy is based on 
the fundamental band for T = 2500°K, pressures up to two atmospheres, and 
turbulent velocities up to 10 km/ sec. The mean absorption coefficient s T m /Av 
has been determined for CO for temperatures of 1000, 2500, and 3500°K. 
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THEORETICAL MOLECULAR LINE ABSORPTION OF CO IN 
LATE SPECTRAL TYPE ATMOSPHERES* 

by 

Virgil G. Kunde 
Goddard Space Flight Center 


INTRODUCTION 

The physical state of early spectral type atmospheres is fairly well known from both an ob- 
servational and a theoretical viewpoint. In general, the opacity sources for these atmospheres are 
considered known, and representative model atmospheres have been developed to the point where 
the effects of non- greyness and line absorption have been included in the models. The current 
central problem in early type model atmospheres is the inclusion of non-LTE effects (Mihalas, 
1966). 

In contrast to the early spectral type atmospheres, the physical state of late spectral type at- 
mospheres is relatively unknown observationally and theoretically. Until recently, little work had 
been done on late type atmospheres mainly because the sources of opacity, especially the molecular- 
line opacity, and their abundances were not known. Recent observations have identified some 
dominant sources of molecular opacity, essentially CO and H 2 0, so that it is now feasible to con- 
struct exploratory late type model atmospheres including the effects of molecular absorption. In 
the construction of late type models, the main problem is knowing how to include molecular opacity 
in the algorithms already developed for non-grey models with continuous opacity sources. Several 
approximate methods for attacking the line-opacity problem have been discussed in the recent liter- 
ature. They are (1) WT V -HT statistical procedure of Strom and Kurucz (1966), (2) the "rectangular- 
box" procedure of Gingerich (1966), (3) the harmonic-mean procedure of Auman (1967), and (4) the 
smeared-line model of Tsuji (1966a, 1966b). Besides line opacity, the equation of state and convec- 
tion also presents problems in constructing realistic late type model atmospheres (Vardya and 
Kandel, 1967: Gingerich, Latham, Linsky and Kumar, 1967). 

The main object of this paper is to find the simplest technique for adequately representing the 
molecular absorption of CO at elevated temperatures. The investigation consists essentially of 
three parts: 

(1) Using an anharmonic vibrating-rotator molecular model, the line positions and line intensi- 
ties have been calculated for the fundamental, first- overtone and second- overtone vibration-rotation 


*Thesis submitted to the Faculty of the Graduate School of the University of Maryland in partial fulfillment of the requirements for the 
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bands of C 12 0 16 and C l3 0 16 . The line intensities have been computed for temperatures from 175 
to 3500°K; collisional-line half -widths have also been determined for each line. 

(2) The CO straight mean mass absorption coefficient has been calculated for several late type 
stellar atmosphere conditions. The mean absorption coefficient has been evaluated by line-by-line 
integration techniques, an approximation S r m /Av, and by the quasi-random band model. The line-by- 
line procedure forms the basis for accurate spectrum analysis and for determining the molecular- 
absorption parameters required by the various methods proposed for solving the line-opacity 
problem, and also checks the accuracy of these same methods. 

(3) The straight mean mass absorption coefficient s T m /Av is determined for the vibration- 
rotation bands of CO at 1000, 2500, and 3500°K. 

A secondary purpose of this paper is to develop the quasi-random band model absorption coef- 
ficient to represent molecular line absorption. This model removes some of the limitations of the 
methods mentioned previously for solving the line opacity problem. The results of this investiga- 
tion show that CO absorption can be represented accurately by the simpler S T m /Lv approximation. 
Thus the quasi-random model has no application in this investigation. 


GENERAL SURVEY OF MOLECULAR ABSORPTION IN LATE 
SPECTRAL TYPE ATMOSPHERES 

Dissociation Equilibria 

It has been known for some time that molecules occur in the solar atmosphere and in the at- 
mospheres of stars whose spectral type is later than the sun's. It is possible to predict the mo- 
lecular abundances in these atmospheres as a function of temperature, gas pressure, electron 
pressure, and initial atomic abundances by assuming the various constituents to be in dissociation 
equilibria. Dissociation equilibria studies (Russell, 1934; de Jager and Neven, 1956; Tsuji, 1964; 
Dolan, 1965; and Vardya, 1966) explain fairly well the observed spectra of the late type stars and 
the branching of the spectral sequence into the K-M, the K-S, and the G-R-N sequences. These 
studies have indicated that the resultant equilibrium concentrations of certain atomic and molecu- 
lar species are extremely sensitive to the initially assumed H:C:N:0 ratio. In general, diatomic 
molecules are predominant in the temperature range 2500 to 5000°K and— with the exception of H 2 , 
CO and N 2 — polyatomic molecules are predominant for temperatures less than 2500°K. The mol- 
ecules which were predicted to be most abundant and also have active vibration-rotation bands are 
listed in Table 1 in order of decreasing abundance for the temperature range of 1000 to 2500^. 
At these low temperatures the outgoing atmospheric radiation is concentrated in the red and in- 
frared. Thus the molecular vibration- rotation bands are an important source of opacity, as they 
absorb strongly in the spectral regions where the outgoing radiation is a maximum. 
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Infrared Observations 


Table 1 


Historically, infrared observations of stel- 
lar objects have been limited by the small 

Molecules in Late Spectral Type Atmospheres 
(Tsuji, 1964). 

amount of infrared energy radiated and by ab- 
sorption in the earth's atmosphere. Recent im- 
provement in infrared-detector sensitivity has 
led to an increasing number of infrared photo- 
metric and spectrophotometric measurements 
through atmospheric "windows." 

Type 

Dwarf 

Giant 

Supergiant 

M-type stars: 

h 2 o 

CO 

ch 4 

nh 3 

h 2 o 

CO 

co 2 

h 2 o 

CO 

co 2 

Several groups have been active in infrared 
photometry. Wildey and Murray (1964a, 1964b) 


CO 2 
HCN 



made the first photometric stellar observations 
through the atmospheric "window" between 8 
and 14 m of 25 stars of spectral type B8 to M7. 
Wildey (1966) extended this work to 42 additional 
stars. A program of multicolor observations 
including 11 broad spectral intervals— of which 

Carbon stars: 

ch 4 

HCN 

CO 

c 2 h 2 

ch 3 

HCN 

ch 4 

CO 

c 2 h 2 

c 2 h 4 

HCN 

CO 

c 2 h 2 

ch 4 

8 occur in the infrared— is being carried out at 






the Lunar and Planetary Laboratory, University of Arizona (Johnson, 1962; Johnson and Mitchell, 
1963; Low and Johnson, 1964; Johnson, Mendoza and Wisniewski, 1965; Johnson, 1966a). These ob- 
servations have contributed to interpretation of the problem of interstellar extinction and have been 
used to revise bolometric corrections and effective temperatures for various spectral types (John- 
son, 1966b). An infrared sky survey is being carried out by Neugebauer, Martz, and Leighton (1965) 
of the California Institute of Technology for the spectral regions 0.68-0.92m and 2. 01-2. 41m. Results 
of the infrared photometric observations have indicated a number of very cool stars, some with 
temperatures apparently ranging as low as 700 to 1000°K (Johnson, Low, and Steinmetz, 1965; 
Johnson, Mendoza and Wisniewski, 1965; Johnson, 1966a). 

Ground-based infrared spectrophotometric observations by Boyce and Sinton (1963, 1964a, 
1964b), Sinton (1966), and Kuiper (1963, 1964) have led to the identification of H 2 0 and CO in the 
atmospheres of M giants and supergiants. On the Stratoscope II balloon flight of November 1963, 
the spectra of six red giants were obtained in the 1 -3 m range (Woolf, Schwarz schild, and Rose, 
1964). The spectra indicated absorption bands of H 2 0 at 1.13, 1.4, 1.9, and 2. 7m and of CO at 2.35m. 
For Mira, it was estimated that H 2 0 absorbs about one-quarter of the total emitted energy. 

Spinrad and Newburn (1965) have obtained low-dispersion spectra of 36 late type stars. The 

o 

(201) vibration-rotation band of H 2 0 at 9400 A was used to identify H 2 0 in these spectra, absorp- 
tion being strong only in the coolest types, M7 and M8. The total amount of H 2 0 for the atmosphere 
of Mira was estimated as approximately 1/3 gm/cm 2 near minimum light. The total H 2 0 abundance 
for Mira was found to vary with phase— the H 2 0 content being reduced as maximum light is ap- 
proached (Spinrad, Pyper, Newburn and Younkin, 1966). Sinton (1966), also, has observed a variation 
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in band strengths with changing phase for the Mira type variables. For R Leo, the intensity of the 
H 2 0 bands were strongest at minimum light and varied with phase appreciably while the CO-band 
intensity was relatively constant. In y Cyg the phase variations of the H 2 0 and CO bands were found 
to vary in opposite directions. For 10 cool stars, Spinrad (1966) has identified several quadrupole 
lines of the first -overtone band of H 2 in the 1-micron region. For Mira and aOri, the data indicated 
a phase change in H 2 abundance. For Mira, the H 2 abundance varied from 6 x 10 26 to 2 x 10 26 
molecules/cm 2 as the phase varied from minimum to maximum light. Spinrad and Vardya (1966) 
have compiled a table of approximate stellar abundances for H 2 0, H 2 , and CO for seven M stars. 

The range in abundances is as follows: 

H 2 0: 2 x 10 22 - 1 x 10 21 molecules/cm 2 , 

H 2 : 6 x 10 25 - 1 x 10 26 molecules/cm 2 , 

CO: 1 x 10 24 - 2 x 10 22 molecules/ cm 2 . 

The above abundances are quite uncertain because of difficulties in observation and data reduction. 

Both the first (Av = 2) and second (Av = 3) overtone bands of CO appear in spectra of M, N, and 
S stars obtained by McCammon, Munch, and Neugebauer (1967) in the spectral regions 1.5-1.8/x and 
1.9-2. 5 fx. The H 2 0 absorption was found to be stronger in the M-type stars than in the N and S type. 

The above infrared photometric and spectroscopic observations qualitatively verify the pre- 
dictions of dissociation equilibria studies that H a O, H 2 , and CO are major molecular constituents 
in late type atmospheres. Considerable quantitative disagreement remains between abundances 
predicted by model-atmosphere theory and observational abundances (Spinrad and Newburn, 1965; 
Spinrad, 1966). Presumably these differences will disappear when more realistic model atmos- 
pheres are constructed which include the effects of molecular absorption and convection. 

Model Atmospheres-Continuous Opacity Sources 

The only set of late type models now in existence has been computed by Gingerich and Kumar 
(1963); Kumar (1964); and Gingerich, Latham, Linsky, and Kumar (1967). These radiative equilib- 
rium models include only the effects of continuous sources of opacity: H, H", H 2 ", He’, electron, 
and Rayleigh scattering from atomic and molecular hydrogen. Molecular absorption and convec- 
tion have not been included in these models. They were constructed for a temperature range of 
2500 to 4000^ and a surface -gravity range of log g = 2-5. They predict several interesting prop- 
erties for late type atmospheres (Gingerich, Latham, Linsky, and Kumar, 1967): 

(1) A large fraction of the hydrogen present is in molecular form. For example, in the outer 
layers of the 2500°K model, over 99 percent of the hydrogen is in molecular form. 

(2) Very high pressures occur in the dwarf photospheres (~200 atm for log g = 5). Inclusion 
of molecular absorption is expected to lower these pressures by 1 to 2 orders of magni- 
tude. Thus, even for these lowered pressures, the pressure-dependence of the molecular 
absorption may have to be included. 
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(3) Convective instability occurs at shallow optical depths, r ~ 0.1. 

Properties (1) and (3) are not expected to change drastically with the construction of more realistic 
models. 


Model Atmospheres-Molecular Absorption 

Only a few investigations concerning molecular -absorption coefficients at elevated tempera- 
tures have yet been recorded. Yamashita (1962) has calculated the molecular absorption of (1) 
metallic atoms and the electronic bands of TiO; (2) rotational lines of CO, OH, and H 2 0; and (3) 
vibration- rotation bands of H 2 0 and C0 2 . The vibration- rotation bands of CO were not included. 

The absorption coefficient was averaged over very wide spectral intervals; for example, the 
vibration-rotation bands were averaged over the entire band. 

To aid in the analysis of the Stratoscope II observations, Auman (1966, 1967) has theoretically 
calculated the molecular absorption of H a O. The positions and strengths of the individual lines 
were calculated for the spectral region from 0.8-12. 5 m for temperatures of 1680, 2016, 2520, and 
3360°K. The harmonic-mean opacity was calculated for spectral intervals of 100 cm" 1 , assuming 
the lines to be Doppler -broadened, for rms turbulent velocities of 2, 4, and 8 km/sec. Water 
vapor was found to be the dominant source of opacity for T < 2520°K; the ratio of the opacity with 
H a O to that without H 2 0 being as large as 2.4 x 10 3 at 2.60/x and 2520°K. The H 2 0 opacity has 
been included in a model atmosphere with T = 2520°K and log g = 3.0. The inclusion of the H s O 
opacity decreased the density at a given temperature in the convective region by a factor of 7 and 
lowered the boundary temperature by 9 percent. 

In a recent investigation, Tsuji (1966a, 1966b) has determined the molecular absorption of the 
vibration-rotation and pure-rotation bands of H a O, CO, and OH. These calculations were based on 
a smeared-line model (Golden, 1962; Patch, Shackleford, and Penner, 1963; and Golden, 1967) to 
approximate the molecular absorption. This model assumes that the rotational structure is smeared 
into a pseudo- continuum either by pressure-broadening or by a sufficient number of rotational lines 
to fill in all spectral gaps. Tsuji also assumes that the H 2 0, CO, and OH molecules can be repre- 
sented by a rigid rotator. On the above foundation Tsuji (1966c) has constructed several preliminary 
model atmospheres of solar composition. 


MOLECULAR LINE ABSORPTION COEFFICIENT 

A molecular absorption line can be described by four parameters: line position (^ 0 ), line in- 
tensity (S), relative line shape (b), and line half-width (a)— see Figure 1. The absorption coef- 
ficient is written in terms of these parameters as 

k(v) = s b(v - V 0 , a) , (1) 
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where the relative line shape b(y-y 0 , a) de- 
pends on the type of line-broadening mechanism 
considered. The line intensity is defined by the 
condition 

/•+00 

S = I k(y) 6.v . (2) 

•'—CO 

Line-by-Line Basis 

The total line absorption coefficient at wave 
number v is 

k(v, P,T) = (T) b(v ~ v o. ■ (P. T)) , (3) 



Figure 1 — Relative absorption coefficient in terms of 
mol ecular -absorption -I ine parameters. 


where the summation includes all the individual absorption lines. As the monochromatic absorp- 
tion coefficients are too numerous to handle routinely the absorption coefficient is averaged over 
a finite spectral interval 


k A v ( p - T) 


1 

Ay 


k(v, P, T)dy . 

J Av 


(4) 


The method used to solve Equation 4 is patterned after the numerical technique developed by 
Drayson (1964, 1966). Basic to this procedure are a variable wave-number mesh, for evaluating 
the monochromatic absorption coefficient, and a division of the contribution of each line into a di- 
rect and a wing contribution. The basic spectral interval considered is Ay = 0.1 cm" 1 . This in- 
terval is divided into a variable number of subintervals depending on the positioning of the line 
centers within the interval. The above integral is evaluated numerically, using Legendre- Gauss 
quadrature (Stroud, 1966), as 


where 


k Av (P. T) = ]2~ Bi / Ai - 22 k K ’ P ’ T ) 

i i=i 


V Bi ~ V Aj "Bj +V Aj 

2 x i + 2 ' 


(5) 


( 6 ) 


The j -summation is for the subintervals in Ay; the i-summation represents four-point Legendre- 
Gauss quadrature for each subinterval, and v B and v k are the upper and lower wave- number 
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boundaries of the subinterval. The function k(v. . , P, t) represents the monochromatic absorption 
coefficient at wave number v. . in A v . 

1 J 

The wave-number mesh v. for a 0.1-.cm -1 interval is constructed in accordance with the fol- 
lowing criterion. A single line is divided into subintervals with the subinterval boundaries placed 
at 0.0, 0.005, 0.01, 0.02, . . . , 0.09, 0.10 cm' 1 from the line center. For lines whose centers are 
closer together than 0.1 cm' 1 , the mesh is adjusted so that each line satisfies the above criterion 
until the midpoint between the two line centers is reached. If no lines occur in a 0.1-cm _1 interval, 
no subdivision occurs and only four mesh points are used. The maximum number of mesh points 
that can occur for an interval of 0.1 cm' 1 is 80. 

For each line, the monochromatic absorption coefficient for distances less than and greater 
than 4 cm' 1 from the line center was considered as the direct and wing contribution, respectively. 
The direct contribution at v.. within an 0.1 cm' 1 interval is 


T ) = ( T )b b-*V“i< P . T >] 


< < V B • ' 4 < "o. < y B + 4 ’ v ij ~ V 0. < 4 cm 1 


(7) 


with the relative line shape given by the mixed Doppler -L or entz line shape. The mixed- line shape 
was evaluated by means of a computer subroutine developed by Young (1965). The wing contribu- 
tion is given by 


K.p.t) = 22 Si ( T)b hj -*v (P ’ T) ] 


- 4 £ 


+ 4 


( 8 ) 


For conditions representative of cool star atmospheres, the line shape is determined by collisional 
broadening for distances from the line center greater than 4 cm' 1 . Substituting the line shape given 
by the Benedict modification of the Lorentz line shape (Equation 47) and the collisional half-width 
of Equation 52, in Equation 8 yields 


k wKi’ p ’ T ) = 22 




v . . - v >4 cm 

ij °i 


(9) 


where the half -width term in the denominator has been neglected. Rewriting Equation 9 gives 


k „ O'ij.P. T) 



(10) 
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where 


SM(v, T) 



S. (T) a C.F. 

1 c 0 


v - v 0 >4 cm 1 . 


(ID 


Since the function SM(y, t) does not vary rapidly with wave number, it is calculated at every 1 cm' 1 ; 
interpolation gives intermediate values. The total monochromatic absorption coefficient is obtained 
by adding the direct and wing contribution: 


kK.P. T ) 


k D ( 


V . . 

1 J 



t) . 


( 12 ) 


Combining Equations 3 and 4 gives, for the mean absorption coefficient, 


k A „ (?. T) 



b ( 


V - V 


0 . 

1 



(13) 


If the spectral interval Av is chosen large enough for the entire absorption line to lie within Ay, the 
mean absorption coefficient is 


(P. T) 




(14) 


where the range of i is for all lines in Ay, and S T is the total line intensity in Ay, 


S 


T 



The mean absorption coefficient S T /Ay is only approximate for intervals where wing effects become 
significant— for example, in regions past band heads. 

At elevated temperatures the number of molecular absorption lines that must be considered— 
for a polyatomic molecule in particular— is so great that even present-day computing facilities are 
not adequate for routine line-by-line computation. Auman (1967), for example, considered approx- 
imately 2.3 million H 2 0 lines in his investigation. Thus it is desirable to resort to a different ap- 
proach; this is discussed in the next section. 
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Band Models 


Various band models have been proposed to represent molecular band absorption for a finite 
spectral interval. Two models, in general, have received the most attention. The Elsasser model 
assumes that the spectral lines are of equal intensity and uniformly spaced; while the random or 
statistical model assumes that the line intensities can be represented by a probability distribution 
with the line positions occurring at random. A third model, the "quasi-random" model, has recently 
been developed to correct certain defects of the Elsasser and the random models. For the purposes 
of this investigation, only the random and quasi-random model are of interest. 

Random Model 

The mathematical formulation of the random or statistical model is developed in Appendix A. 
Assuming the line positions occurring in spectral interval At> to be distributed at random, Equation 
A5 gives the average absorption coefficient k(v) at v due to n lines in Av as 


kO) 




• . a. j P' (a.) P(S.) dv 0 . 

P ( S i) p ' K) dl V dS i d “i 


dS. da. 
1 1 


(15) 


Quasi-Random Band Model 

To correct certain limitations of the Elsasser and statistical band models, the quasi-random 
band model has been developed (Wyatt, Stull, and Plass, 1962). The quasi-random model differs 
from other band models in several ways: 

(1) Method of spectral division: To determine the average absorption coefficient for the in- 
terval A, A is subdivided into smaller intervals e. The vibration-rotation lines occurring in e are 
assumed to have a random distribution of line positions. The average absorption coefficients for 
the e's within A are then arithmetically averaged to give the average absorption coefficient for A. 

In principle, the exact line positions are obtained by allowing e to approach zero. 

(2) Line -intensity distribution: The lines in a subinterval e are divided into subgroups, each 
subgroup covering one intensity decade. The intensity distribution is then simulated by the number 
of lines and the average intensity s for each decade. For four intensity decades, Equation 15 yields 
the average coefficient due to one subinterval as 


4 /• 

k(v) = 2^7 J n t Sj b(v - v 0 _ , dv„. , (16) 

1 = 1 e 


II 



where the total number of lines in e is 


4 



1 = 1 


(3) Wing effects: The contribution to the absorption coefficient of a given subinterval from 
the wings of lines in adjacent subintervals is properly taken into account. 

For the quasi-random model, the total average absorption coefficient at v in subinterval e k is 


k k Ol = £>, (v)l (17) 

i = 1 


where k. (v) is the contribution to k, (v) from the n lines in interval e.. The direct contribution 
k k ( v ) is caused by the lines in interval <r k when j = k. The wing contribution is due to all remain- 
ing intervals for which j / k. Combining Equations 16 and 17 gives 





( 18 ) 


If the e interval is chosen large enough for the lines to lie entirely within e, Equation 18 reduces to 
S T /Av, which is the same as the expression given by Equation 14. From Equation 18, the average 
absorption coefficient for each subinterval e can be evaluated. The quasi-random average absorp- 
tion coefficient for interval A is determined by arithmetically averaging the k k 0)'s of the subin- 
tervals contained in A. The line shape is treated in the manner described in the next section. The 
integral in Equation 18 is evaluated numerically. 

To evaluate Equation 18 requires the specification of the following parameters: 

(1) € , A, subinterval and interval size (in this investigation e = 5 cm" 1 , while A has been varied 
from 5 to 100 cm' 1 ). 

(2) nj, Sj, aj— the number of lines, the average line intensity, and the average collisional-line 
half -width for each intensity decade and for each subinterval. 

(3) Line shape, b^-v 0 , ajj. 

A program was developed for the IBM 360, to determine values of nj , Sj , and a c for each 
subinterval, from input tables of the individual line parameters v 0 , S, and . 


s. 
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MOLECULAR LINE PARAMETERS FOR CO 


Line Position 

The term value of a vibrating anharmonic rotator (Herzberg, 1950) is 

T' (cm' 1 ) = (v + 1/2) - x e (v + 1/2) 2 + co e y e (v + 1/2) 3 

+ z e (v + 1/2) 4 + B v J(J + 1) - D y J 2 (J + l) 2 , (19) 


where 


B v = B e - a e (v+ 1/2) + r(v+ 1/2) 2 , (20) 

D v = D e + /3(v + 1/2) . (21) 

It is convenient to define a rotational quantum index m by m - J + l for the R branch and m = - J for 
the P branch. The term values given by Equation 19, plus the vibration and rotation selection rules, 
yield the predicted line positions. The constants needed to solve Equation 19 are well known for 
CO. The values listed by Benedict, Herman, Moore, and Silverman (1962) were used in this in- 
vestigation, with the exception of for C l3 0 16 which is incorrect in their paper. The correct 
value is 2121.462 (Mills and Thompson, 1953). These constants yield line positions that agree with 
experimental line positions within 0.05 cm' 1 and 0.5 cm' 1 for the fundamental and first-overtone 
bands, respectively. 

Line Intensity 

An optical depth can be defined by 


dr - - k m ( v ) ds . (22) 

where k m (v) is the mass absorption coefficient, p s is the density of the absorbing gas, and s is the 
geometrical path length. To conform to the system of units used in laboratory measurements, it is 
necessary to transform the mass absorption coefficient k m (v) into an absorption coefficient per unit 
length per unit pressure k lp (v) . The transformation is made by defining 

k lp (v) = k m ( V )p g /p, (23) 

and 

du(cm atm) = - p ds , (24) 
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where p is the partial pressure of the absorber in atmospheres and u is the optical path length in 
cm atm. The absorption coefficient k lp (v) has units of cm -1 atm -1 . Equation 22 may now be 
written 


<Jt - k lp (y)du . (25) 

The absorption coefficient transformation also requires redefining line intensity. The normaliza- 
tion condition for the line strength requires that 

Z' + OO 

k m O)dv = S m (cm - gm _1 ) , (26) 


and 



(v) &v 



(27) 


where the transformed line intensity S lp is related to S m by 

, S m M(l. 01325 X 10 6 ) 

S lp = sm V p = R T 

m 


(28) 


where R = 8.314 x 10 7 , and M is the molecular weight. The factor 1.01325 x 10 6 is due to the 

m ' 

change in pressure units. 

The line intensity, S lp , for a transition from quantum state 1 to quantum state u is given by 
Penner (1959, Equations 7-124) 


s ip 


877 3 V . 

lu 

3 he gj 


( N i/P) 



1 

Ir 

n 

l R lJ 2 

1 - exp - ^ kX 


(29) 


where v is the transitional wave number, gj is a statistical weight, is the number of molecules 
per unit volume in the lower quantum state, and R lu is the electric-dipole-moment matrix element. 
Equation 29 is based on a distribution of quantum states in thermodynamic equilibrium. The dipole- 
moment matrix elements are related to the Einstein coefficient for induced absorption (Penner, 
1959) by 


B 


lu 



(30) 
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Penner defines the dipole-moment matrix element for a transition from a lower state of vibrational 
and rotational quantum numbers v and J to a upper state of v' and J' as 


K 

1 2 = s T Ir"',- 1, 1 2 . 

u 1 J 1 V ,J 1 ’ 

(31) 

R v v y r = 

[ >A v » i j< ( r ) M (0 '/'v.J ( r ) dr • 
•'0 

(32) 

s j 

- J , P branch , 

(33) 

s j 

= J + 1 , R branch , 

(34) 


where r is the internuclear distance, M( r ) is the dipole moment, and \p( r ) represents the solution to 
the radial Schrodinger equation. Using the Maxwell -Boltzmann distribution and Equation 31 trans- 
forms Equation 29 to 



where N 0 is the total number of molecules per unit volume and Q is the partition function. The 
above expression is normalized for standard temperature (T 0 = 273.18°K) and pressure (p 0 = 1 
atm) by 


o lp 

°STP 


(T) 


N 0 (T 0 ) p 
Po N 0 (T) 


S lp (T) 



(36) 


From Equations 35 and 36, 


Srtp (T) 


877 3 V, 

l_u 

3 he 


N o (Tq) 
Po 



7.J 



(37) 


Assuming the ideal gas law, 



(38) 


This transformation allows the line intensity to be computed as a molecular parameter independent 
of the actual number of molecules present. In this normalized system of units, the optical path 
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IS 


du gTP (cm atm) - (T o y'T)du 


(39) 


and the optical depth is 


A t — If 1 P Hi | 
u/ k stp stp 


(40) 


The corresponding transformation to mass absorption units for the line intensity is 


S m (T) = 


R T n 

m 0 

M (l x 10 6 ) 


Ss'/p (T) • 


(41) 


The dipole moment is expressed as a power-series expansion about the equilibrium inter- 
nuclear separation, r e : 


M(r) = M i ( r “ "e) 1 ‘ 


(42) 


The expansion coefficients M. are determined from laboratory measurements of infrared band in- 
tensities. By this method, Young and Eachus (1966) have determined the dipole moment of CO to 
be 


M(r) = + .112 + 3.11 (r -r e ) - .15 (r- r e ) 2 - 2.36 (r- r e ) 3 . (43) 

The dipole moment matrix elements of Young and Eachus, calculated by numerical evaluation of 
the integral of Equation 32, were used in this investigation. The accuracy of the matrix elements 
determined in this manner depends mainly on the accuracies of the infrared band intensities and 
the Taylor expansion of the dipole moment. 

The partition function is given by Herzberg (1950): 


Q 


LL 

„ = n t = n 


(2J + l)e 


,i/ kT 


and is presented in Figure 2 for C 12 0 16 and C 13 0 16 . 


(44) 
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Line Shape 


A -monochromatic absorption line may be 
broadened by one or more of many different 
types of line -broadening mechanisms, both 
physical and geometrical. The mechanisms of 
importance for this investigation are 

. (1) Collisional or pressure broadening 
caused by collisions with neighboring 
atoms or molecules, 

(2) Doppler broadening caused by random 
molecular thermal motions, 

(3) Turbulent broadening. 

For a collision-broadened Lorentz line, 


b (v- ^0. “ c ) 


a 


c 
7 T 



(45) 


where a is the collisional half-width. The 

c 

relative Doppler line shape is 


b{v -v 0 , a D ) 


yin 2 





Figure 2— Partition function of C 12 O 16 and C 13 O 16 
as function of temperature. 


where a D is the Doppler half -width at half-maximum. As the Lorentz line shape is theoretically 
valid only near the line center (1-3 cm -1 ) (Bohm, 1960; Goody, 1964), various experimental studies 
have been made to determine the line shape of the collision-broadened wing. The experimental 
studies for various molecules such as HCL, CO, and C0 2 have indicated that the line wings are 
not Lorentzian. Experiments on absorption by the fundamental and first overtone of CO (Benedict, 
Herman, Moore and Silverman, 1962) have shown that the Lorentz line shape should be modified 
by an exponential multiplying factor, dependent on \v - v 0 1 for | v - v Q 1 > d cm -1 . The modified line 
shape is 


b(v-^ 0 ,a c ) - b£j° t ' z ac) C.F. , |v~v 0 | > dcm' 1 , (47) 

where the correction factor (C.F. ) is 

C.F. = e" a (l V ’ / °l' d ) b . (48) 

The line shape of Equation 47 is called the "Benedict modification" of the Lorentz line shape. The 
constants were found to be d =4 cm -1 , b' = 1, and a = 0.015 for CO- CO collisions. In this study 
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it is assumed that the CO line shape caused by broadening by foreign gases is the same as for 
self-broadening. 


The line half -width caused by thermal and turbulent broadening was assumed to be 


a 


, /2kT ; 

y/c y nr ln 2 + • 


(49) 


where £ is the rms turbulent velocity. 


Collisional-Line Half-Width 

Kinetic theory relates the line half-width to the effective collision cross-section a (Smith, 
Lackner, and Volkov, 1955) as 


a 


c 


N / .(T) V(T ) cr(T) 
2tt c 


(50) 


where N' (T) is the number of particles per unit volume and V(T) is the mean relative velocity of 
the colliding particles. Theories for determining the effective collision cross section have de- 
veloped along several lines. For optical transitions, collisions may be described adiabatically; 
in the infrared and microwave region, diabatic collisions (in which the collision can induce a 
transition in the radiating molecule) must be considered. The diabatic theory of Anderson (1949) 
allows the effective collision cross section to be defined in terms of the intermolecular forces 


This theory is most easily applied when one type 



Figure 3 — Collisional half-width of CO as function of rotational 
quantum index at pressure of one atmosphere and room tempera- 
ture. The experimental values are for CO-H 2 collisions (Crane- 
Robinson and Thompson, 1962b). The solid curve was adopted for 
this investigation. 


of intermolecular force dominates. The largest 
collision cross sections are associated 
with long-range forces (dipole-dipole). 

In general, experimental cross sections 
are larger than would be estimated from 
kinetic theory (Townes and Schawlow, 
1955). 

Collisional line widths have been 
measured for self-broadenedand foreign- 
broadened CO (Eaton and Thompson, 1959a, 
1959b; Crane -Robinson and Thompson, 
1962a, 1962b; Benedict, Herman, Moore 
and Silverman, 1962). Figures 1 and 2 
of a paper by Shaw and Houghton (1964) 
plot the experimental-line collisional 
half-widths of the fundamental band of 
C 12 O ie for CO-CO and CO-N 2 collisions 
as a function of rotational quantum num- 
ber. Figure 3 of the present paper (Kunde) 
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plots the experimental half-widths for the CO fundamental for CO-H 2 collisions for T = 300°K and 
P = 1 atm. The experimental data are available only for low values of the rotational quantum num- 
ber. For application to cool star atmospheres, the half -widths must be known for higher values of 
the rotational quantum number; also, the half-width dependence on temperature must be known. 

Figure 3 also presents the line half-widths corresponding to a kinetic cross section, to London 
dispersion forces, and to a hard-sphere collision. The observed cross section is evidently smaller 
than the kinetic cross section. It has been suggested that this is because of short-range interac- 
tions that depend on a high power of 1/r (Anderson, 1949). Long-range interactions are not ex- 
pected to be important in CO-H 2 collisions, as H 2 has no dipole moment and CO has only a very 
weak dipole moment. For CO-H 2 collisions, the short-range interactions that may possibly be 
important are 

(1) dipole -quadrupole ( a r -4 ), 

(2) quadrupole -quadrupole (<* r“ s ), 

(3) dipole -induced dipole ( a r' 6 ), 

(4) London dispersion (ar -6 ), 

(5) quadrupole -induced dipole (« r -7 ), 

(6) exchange forces due to direct interaction of the electronic-charge distributions of the col- 
liding molecules. 

Crane-Robinson and Thompson (1962b) have evaluated types (3) and (4) and find London dispersion 
forces the most significant. In Figure 3, the half -widths due to London dispersion forces do not 
have the observed dependence on rotational quantum number. The half-widths given by the solid 
line were used where the half-width for high J values was somewhat arbitrarily assumed to cor- 

o 

respond to a hard-sphere collision diameter of o- = 2 A. This same curve is assumed to apply to 
all vibrational levels. Penner (1952) has shown that the cross section is not a sensitive function 
of vibrational level. 

The temperature dependence of the line half -width for interactions of the type l/r n is 


a « 1 



(51) 


if n is large (Townes and Schawlow, 1955). From Equations 50 and 51, the collisional-line half- 
width for CO-H 2 collisions at pressure Pand temperature T can now be written as 


a c„ (P 0 . V m ) P/Po/yr ’ 


(52) 


a c (P, T, m) 
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where a c (p 0 , T 0 , m) is given by the solid curve of Figure 3. In any event, the accuracy of 
the collisional half-width is not a critical factor in the stellar-atmosphere situations con- 
sidered in this investigation; the Doppler and turbulence broadening are the dominant line- 
broadening mechanisms. 


Using the formulation and basic molecular data discussed above, a computer program was 
developed for an IBM 360 to generate line position, line intensity, and collisional half -width at a 
given temperature for each line of the fundamental and first- and second- overtone bands con- 
sidered. The C 12 0 16 andC 13 0 16 line intensities were weighted according to a terrestrial C 12 O 16 / 
C 13 0 16 abundance ratio of 89.2 (Benedict, Herman, Moore and Silverman, 1962). Tables of line 
positions, line intensities, and collisional half-widths are given for a temperature range from 175 
to 3500°Kfor the fundamental and first- and second-overtone bands (Kunde, 1967). Figures 4, 5, 
and 6 present the line intensities, S m , at 2500°K for the fundamental and first- and second- overtone 
bands, respectively. For a given band at a given temperature, the calculations include all lines 
within approximately 5 to 6 orders of magnitude of the maximum line intensity. At 2500 K, 3030 
lines were considered for the fundamental band. Young (1966) has presented the line intensities, 
Sg'/p, for the fundamental and first- and second- overtone bands at T = 3000°K. 




Figure 4 — Line intensities corresponding to a mass ab- 
sorption coefficient for various vibrational levels of the 
C ,2 0 16 fundamental band at 2500°K. 


Figure 5 — Line intensities corresponding to a mass ab- 
sorption coefficient for various vibrational levels of the 
C 12 O 16 first -overtone band at 2500°K. 
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COMPARISON OF THEORETICAL AND 
EXPERIMENTAL ABSORPTION DATA FOR CO 

To check the accuracy of the CO molecular 
parameters computed in the previous section, 
the absorption by CO has been computed the- 
oretically by direct integration. The theoreti- 
cal absorption was then compared to existing 
laboratory absorption data at planetary tem- 
peratures (Burch, Gryvnak, Singleton, France, 
and Williams, 1962) and at elevated tempera- 
tures (Abu-Roma and Tien, 1966). Figures 7 
and 8 present comparisons at 300°K for the 
fundamental and first-overtone bands, respec- 
tively. The figures also indicate the experi- 
mental pressures, optical path lengths, total 
band absorptions (W) , and the theoretical total 
band absorptions. The theoretical computations 
are based on the Benedict modification of the 
Lorentz line shape. The theoretical and ex- 
perimental absorption curves are indicated by 
the dashed and solid curves, respectively. The 
theoretical transmittances are averaged over 
25-cm" 1 intervals, while the experimental 
transmittances correspond to a resolution of approximately 25 cm" 1 and 15 to 20 cm" 1 for the 
fundamental and first-overtone bands, respectively. The theoretical absorption curves agree with 
the experimental absorption curves almost within the experimental error, although there are in- 
dications that the theoretical curves systematically underestimate absorption in the center of the 
P and R branches. The total band absorptions agree within 10 percent. 

Abu-Romia and Tien (1966) have measured the spectral absorption and the total band absorp- 
tion of the fundamental and first-overtone bands of CO at temperatures up to 1500°K. Table 2 com- 
pares the theoretical and experimental total band absorptions for the fundamental band at 1500°K. 
There is general agreement to within 15 percent except for the first case which can probably be 
attributed to experimental inaccuracy. Comparison of the spectral absorption curves indicates 
that the theoretical absorption is generally greater than the experimental absorption by about 10 
to 20 percent for the spectral range 1950 to 2225 cm" 1 . Considerable disagreement exists for 
wave numbers less than 1950 cm" 1 and greater than 2225 cm' 1 . Possibly, C0 2 absorption affects 
the experimental values on the high wave-number side of the band. 

Figure 9 shows the band intensities as a function of temperature. The solid line results from 
a summation of the line intensities of the individual vibration- rotation lines computed in this in- 
vestigation. The dashed curve represents the theoretical temperature dependence computed for a 



Figure 6 — Line intensities corresponding to a mass ab- 
sorption coefficient for various vibrational levels of the 
£12 o' 6 second -overtone band at 2500°K. 
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1950 2000 2050 2100 2150 2200 2250 

WAVE NUMBER (cm- 1 ) 



WAVE NUMBER (cnf 1 ) 

Figure 7— Comparison of theoretical absorption curves and total band absorption with experimental data 
(Burch, Gryvnak, Singleton, France, and Williams, 1962) for the fundamental band. 


-harmonic-oscillator/rigid-rotator model by Breeze, Ferriso, Ludwig, and Malkmus (1965). 

The solid and dashed curves are normalized at 273 °K. Deviations of the solid and dashed 
curves at the higher temperatures probably result from anharmonicity and vibration-rotation 
interaction. Breeze and Ferriso (1965) have shown that the relative experimental temperature 
dependence of the fundamental and first -overtone CO bands also agrees with the relative the- 
oretical temperature dependence for the range 2500 to 5000°K. Thus, indirectly, the computed 
band intensities of this investigation have the correct relative temperature dependence. 

Comparison of the theoretical and experimental absorption data for CO indicates that the CO 
molecular parameters computed in this investigation predict CO absorption within 10 to 20 per- 
cent at the elevated temperatures. Improvement of molecular-parameter accuracy depends on 
more complete knowledge of the dipole moment for high values of internuclear separation and of 
the collisional line shape for foreign-broadened gases. In addition, more accurate experimental 
absorption data are needed at the elevated temperatures. 
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Figure 8— Comparison of theoretical absorption curves and total band absorption with experimental 
data (Burch, Gryvnak, Singleton, France, and Williams, 1962) for the first-overtone band. 


Table 2 

Comparison of Theoretical and Experimental (Abu-Romia and Tien, 1966) 
Total Band Absorption for CO Fundamental Band at 1500° K. 


u . e q 0 (cm atm) u,„(cmatm) P (atm) (cm 1 ) W,. (cm x ) 




COMPARISON OF EXACT, APPROXIMATE S T m /Av, 
AND QUASI-RANDOM BAND MODEL MEAN 
ABSORPTION COEFFICIENTS 

Exact Mean Absorption Coefficient 

The straight mean absorption coefficient 
was evaluated by direct integration for spectral 
intervals of 0.1, 1.0, and 5 cm' 1 , as shown in 
Figures 10, 11, and 12 respectively. Three 
cases have been chosen to represent late-type 
atmospheres: 

(1) T = 2500°K, P = 0.0 atm, £ t = 0.0 km/ 
sec, Doppler 

(2) T = 2500TC, P = 2.0 atm, = 0.0 km/ 
sec, Collisional-Doppler 

(3) T = 2500°K, P = 0.0 atm, <f t = 10 km/ 
sec, Doppler -Turbulent 

For case 2, the Benedict modification of the 
Lorentz line shape was used. Test computa- 
tions with the Lorentz line shape agreed to 
within 1 percent of the Benedict modification in 
the band center where the wing contribution is 
very small. In regions where the wing contri- 
bution is an appreciable fraction of the total 
absorption coefficient, such as the region short 
of 1700 cm" 1 and the region past the 0-1 band 
head, the Lorentz contribution may be 10 to 100 
times larger than the Benedict modification. 


300 


250 


FUNDAMENTAL 


SUMMATION OF LINE INTENSITIES 


THEORETICAL 

_i i I i i__i l-i i 




Figure 9— Comparison of band intensities as function of 
temperature. The solid lines were obtained by summa- 
tion of the individual vibration-rotation line intensities 
computed in this investigation. The dashed lines repre- 
sent the theoretical temperature dependence computed 
for a harmonic-oscillator/rigid-rotation model (Breeze, 
Ferriso, Ludwig, and Malkmus, 1965). 


The present computation scheme establishes approximately 175,000 wave-number mesh points 
within the fundamental CO band for evaluation of monochromatic-absorption coefficient. The run- 
ning time for one case is approximately 90 to 120 minutes on an IBM 360. The present direct- 
integration program can be optimized considerably for computations at conditions typical of late- 
type atmospheres; say, by approximating the mixed-line shape (rather than using a general 
subroutine) and reducing the number of subintervals for a given line. These modifications would 
probably more than halve the computing time. 

Figure 10 plots the mean absorption coefficient, for = 0.1 cm' 1 , at every 0.1 cm' 1 for the 
region near the center of the R branch. At this resolution, the increasing wing contribution is evi- 
dent (as turbulent or collisional broadening is considered), with the absorption coefficient varying 
by orders of magnitude between the various cases. Figure 11 plots the mean absorption coefficient, 
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Figure 10— Straight mean mass absorption coefficient for 
0.1-cm" 1 intervals for a portion of the CO fundamental 
band at 2500°K. The mean absorption coefficient was 
obtained by direct integration. 


for A v = 1.0 cm" 1 , at every 0.5 cm” 1 . Only the 
cases for pure Doppler and pure collisional 
broadening are indicated. Absorption- coefficient 
variations due to turbulent broadening are of 
the same order as the variations due to pure 
collisional broadening. Figure 12 plots the mean 
absorption coefficient, for Av = 5 cm” 1 , at 
every 5 cm” 1 ; note that the ahsorpti on coeffi- 
cient is on a linear scale. Only case 1 (pure 
Doppler broadening) is shown, as the absorp- 
tion coefficient for the other two cases agrees 
with this to within 5 percent. It is evident that 
as the spectral interval width is increased 
(Figures 10 to 12) the mean absorption coef- 
ficient approaches the limiting condition of 
S^/Av. Thus pressure effects connected with 
the collisional line shape may be neglected for 



Figure 11— Straight mean mass absorption coefficient for 
1.0-cm” 1 intervals for a portion of the CO funda- 
mental band at 2500°K. The mean absorption coefficient 
was obtained by direct integration. 



WAVE NUMBER (cm” 1 ) 


Figure 12 — Straight mean mass absorption coefficient for 
5-cm” 1 intervals for a portion of the CO fundamental 
band at 2500°K. The mean absorption coefficient was 
obtained by direct integration. 
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CO even for pressures of several atmospheres. The negligible pressure effect can be attributed 
mainly to the small CO-H 2 collision cross section. For larger collision cross section, pressure 
effects' would presumably have more effect on the mean absorption coefficient at resolutions com- 
parable to Aw = 5 cm' 1 . Turbulent velocity effects, also, are negligible for Aw > 5 cm ' 1 . 


Quasi-Random Band Model Mean Absorption Coefficient 


The smallest subinterval spectral width chosen for the quasi-random model was e = 5 cm' 1 . 
Comparison of the quasi-random mean absorption coefficient with the mean obtained from direct 
integration for the CO fundamental indicates that agreement is generally within 5 percent, with 
occasional disagreement up to 10 percent for each of the three cases. At Aw = 10 cm' 1 , the agree- 
ment is generally within 1 percent, with occasional disagreement up to 5 to 10 percent. The com- 
puting time for one case, using the quasi-random model, is 1 to 2 minutes on the IBM 360. As il- 
lustrated in the previous section for typical late type stellar atmosphere conditions, wing effects 
are negligible even for Aw = 5 cm' 1 . Under these conditions the simpler S T m /Aw approximation can 

be used and thus the quasi-random model will be ,„4 

10 I ' r ~ r T ' I 1 I 1 I 1 ! 1 I 1 I ' ! 1 I ' 1 1 I 1 I 1 ! ' J 

discussed no further in this investigation. - H 


Approximate S T m /Av Mean 
Absorption Coefficient 

The approximate mean absorption coeffi- 
cient Sy /Aw has been determined for the funda- 
mental band of CO at temperatures of 1000, 
2500, and 3500°K and is shown in Figure 13 for 
Aw = 100 cm' 1 . The dashed portions of the 
curves indicate the region beyond the 0-1 band 
head, where S T ra /Av is only approximate. Com- 
parison of S T m /Aw with the mean absorption co- 
efficient determined by direct integration in- 
dicates agreement with each of the three cases 
within 1 percent, except for the region short of 
1600 cm' 1 where differences of 5 to 10 percent 
occur. Because line wing effects are negligible, 
the dashed portions of the curves in Figure 13 
do accurately represent the mean absorption 
coefficient. Figures 14 and 15 show values of 
S T m /Aw for 100 cm' 1 for the first- and second- 
overtone bands at 1000, 2500, and 3500°K. Even 
at this low resolution, band heads are apparent 
in both overtones. 



Figure 13— Straight mean mass absorption coefficient for 
100-cm" 1 intervals for the CO fundamental band at 
1000, 2500, and 3500°K. The mean absorption coef- 
ficient is based on the S^"/Aw approximation. The 
dashed portions of the curves indicate the region beyond 
0 - 1 band head, where S T m /Aw is only approximate. 
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Figure 14 — Straight mean mass absorption coefficient for 
100-cm" 1 intervals for the CO first-overtone band at 
1000, 2500, and 3500°K. The mean absorption coef- 
ficient is based on the S T m /Av approximation. The 
dashed portions of the curves indicate the region beyond 
0 — 1 band head, where S T m /Av is only approximate. 

CONCLUSIONS 



Figure 15— Straight mean mass absorption coefficient for 
100-cm -1 intervals for the CO second -overtone band at 
1000, 2500, and 3500°K. The mean absorption coef- 
ficient is based on the S T m /A v approximation. The 
dashed portions of the curves indicate the region beyond 
0 — 1 band head, where S T m /Av is only approximate. 


The main results and conclusions of this investigation are: 

(1) From an anharmonic vibrating-rotator molecular model, the molecular absorption line 
parameters have been derived for C 12 0 16 and C 13 O ie . The line positions, line intensities, and 
collisional-line half -widths have been computed for temperatures from 175 to 3500°K. 


(2) The CO straight mean mass absorption coefficient has been evaluated, by direct integration 
techniques, by the approximation S T m /a v , and by the quasi-random band model. The direct integration 
technique allows a check on the accuracy of the S T m /Av approximation and of the quasi-random band 
model. 

(3) The quasi-random model is more accurate than the S T m /Av approximation when the effects 
of line wings must be considered. For late type atmospheres, wing effects were not important for 
pressures up to two atmospheres or turbulent velocities up to 10 km/ sec. The S T m /Av approxima- 
tion represents CO absorption very accurately in late type atmospheres for spectral intervals 
larger than 5 cm ~ 1 . 
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(4) The parameter s T m /hv has been determined for the CO vibration- rotation bands for tem- 
peratures of 1000, 2500, and 3500°K at spectral intervals of 100 cm" 1 . 

Further investigations will extend similar calculations to other molecular absorption bands. 
Above 3000°K, electronic spectra should be important; in the range of 1000°K for M-type stars, 
polyatomic molecules such as CH 4 , NH 3 , and HCN should be considered. For carbon-rich stars, 
the molecular absorption of CH 4 and HCN is important. 
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Appendix A 


Statistical Band Model 


The statistical model (Goody, 1952; Plass, 1958) considers a spectral interval of width t\v con- 
taining n total lines with each individual line having intensity S half -width a. and center wave 
number v. occurring within Av. Let N (v 1 ••• v n ) dv 1 ■ • ■ dv n represent the probability that the i th 
line occurs between- v. and v i + dv. , P(S i ) dS. the probability that the i th line has an intensity be- 
tween s. and s. +dS. , and P' (a. • ■ ■ <x n ) the probability that the i th line has a half-width between 
a. and O.. + da . . Assuming that v 1 ■ ■ ■ v n , • a n , S. • • • S n are mutually independent, the prob- 
ability of finding the set of n lines with ■ • • v n , a 1 ■ ■ ■ a nJ Sj • • • S n is 

N K ••• v n ) P ( S i ••• S n) P ' ( a l a n) dv X ••• dv n dSj ••• dS n da t ••• da n . (Al) 


The average absorption coefficient k(v) at v is found by averaging k(v) over the probability dis- 
tribution of the set of n lines 


k(v) 


-CO -CO -CO -CO - /- 

... ... ... k(v)N(v 1 -"v n )p(S 1 "-S n )P'(a 1 ---a n )d^ 1 ---dv n dS 1 ---dS n da 1 ---da n 

•^0 ->0 -Jo Jo JAv JAv . 

-CO /-CO -CD -CD — — 

... ... ••• N(v 1 ---v)P(S 1 -"S ll )P'(a 1 ---a ii )dv 1 -“dv n dS 1 -"dS n 

Jo Jo Jo Jo Ja* 


(A2) 


da , • ’ ’ da 
l n 


The absorption coefficient k(>) is given by 


k(v) = Z] Si h ( v ~ v ^ ’ 

i = i 


(A3) 


where S. is the line intensity subject to the normalization condition 


— +co 

k i 

■/—CO 


( v ) dv - S . . 


(A4) 


31 


and b (v - v. , a.) is the relative line shape. Assuming that the line positions are distributed at ran- 
dom in Av (v ; ) = constant), k(v) can be written: 


k(^) 



, a.) P' (a.) P(S.) dv. dS. da. 
(a.) P(S.) dv. dS. da. 


(A5) 


For example, consider the case where all the spectral lines are of equal intensity S and equal half- 
width a. Then 


kO) 



nk'W, 


(A6) 


where 


k' (v) 



(A7) 


The parameter k' (v) represents the average absorption coefficient at v— averaged with respect to 
center wave number— for a line of intensity S and half-width a whose center wave number v. occurs 
at random in Av. A second example is given by 


P(S) 



(A8) 


P' (a) - constant , 


(A9) 


for which 


k(v) = 


n S„ 



Av 


a) dv. 


(A10) 
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